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New constraints on the Polarization of Anomalous 
Microwave Emission in nearby molecular clouds 
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ABSTRACT 

Anomalous Microwave Emission (AME) has been previously studied in two well-known 
molecular clouds and is thought to be due to electric dipole radiation from small spin- 
ning dust grains. It is important to measure the polarization properties of this radiation 
both for component separation in future cosmic microwave background experiments 
and also to constrain dust models. We have searched for linearly polarized radio emis- 
sion associated with the p Ophiuchi and Perseus molecular clouds using WMAP 7- year 
data. We found no significant polarization within an aperture of 2° diameter. The up- 
per limits on the fractional polarization of spinning dust in the p Ophiuchi cloud are 
1.7 %, 1.6 % and 2.6 % (at 95 % confidence level) at K-, Ka- and Q-bands, respectively. 
In the Perseus cloud we derived upper limits of 1.4 %, 1.9 % and 4.7 %, at K-, Ka- and 
Q-bands, respectively; these are similar to those found by Lopez-Caraballo et al. If 
AME at high Galactic latitudes has a similarly low level of polarization, this will sim- 
plify component separation for CMB polarization measurements. We can also rule out 
single domain magnetic dipole radiation as the dominant emission mechanism for the 
20-40 GHz. The polarization levels are consistent with spinning dust models. 

Key words: radiation mechanisms: general - polarization - radio continuum: ISM: 
clouds - ISM: individual objects: p Ophiuchus molecular cloud; Perseus molecular 
cloud - cosmology: diffuse radiation 



1 INTRODUCTION 

Anomalous Microwave Emission (AME) has been ob- 
served by numerous experiments in the frequency range 
~ 10-100 GHz as an excess compared to synchrotron, 
free- free, cosmic mi crowave background (CMB) and ther - 
mal dust emissions (ILeitch et al. I ll997l : iKogut et all 1 19961 ) . 
The A ME is observed as diffuse e mission at 
itudes (Ide Oliveira-Costa et all |2004|: 
Miville-Deschenes et all 120081 
20ld: iGold et al 



Ysard et all 



Bonaldi et al 
Dobl er fc Finkbeine 
20111) as well 




spe- 



cific Gal^£H£_obJects ([F^^eine^e^^||200^_^^rtso^^^^] 
20051 ; ICasassus etlll bood 120081: IDickinson et all l2009a| . 



AMI Consortium et all l2009f ; iTibbs et al 
Castellanos et all l201ll : IVidal et all 1201 ll ). The 



201C : 
physical 

mechanism responsible for the AME is a matter of debate. 
However, there is significant evidence fo r electric dipole ra- 
diation from small spinning dust grains dDraine fc Lazarianl 
Il998allb|). In particular, the rec ent results from Planck 
l|Planck Collaboration et alll201ll ) have shown striking evi- 
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dence for a spinning dust spectrum for both the Perseus and 
p Ophiuchi molecular clouds. 

The exact amplitude of the AME signal is still not well 
understood since the results depend critically on the com- 
ponent separation of multiple diffuse components (free-free, 
synchrotron, thermal dust and CMB). For example, in the 
case of template fitting, where AME is correlated with far- 
infrared (FIR) maps, the AME is the strong est foreground 
in total-intensity in the range ~ 20-40 GHz ([Bandav et all 
120031 : iDavies et al]|200d IDickinson et al.ll2009bl ). It is there- 
fore crucial to understand this component in terms of CMB 
foreground removal. Given that the diffuse polarized fore- 
grou nds are much brighter than the CMB polarization sig- 
nal l|Page et all 120071 : iBetoule et al J 120091 ) ■ it is critical to 
understand the nature of the AME and determine its polar- 
ization properties. Although we expect synchrotron radia- 
tion to be the dominant polarized foreground at frequencies 
< 100 GHz, polarized AME could be significant. 

To date, there have been very few measurements of 
the polarization of AME. The first measurements were 
made with the COSMOSOM AS experiment at 11 G Hz in 
the Perseus Molecular cloud. Battistelli ct al. (2006) mea- 
sured the fractional polarization to be 3.4l 19 % at 95 % 
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Figure 1. WMAP 7-year maps of the p Ophiuchi region at 22.7 (K-band), 33.0 (Ka-band), 40.7 (Q-band), 60.6 (V-band), 93.4 GHz 
(W-band). Each map covers 10° X 10° centred at (I, b) = (353.°05, +16.°90) and is smoothed to 1° resolution. From left to right are Stokes 
I (total-intensity), Q, U, polarized intensity (P) and polarization angle (PA). Units are thermodynamic (CMB) mK. The graticule has 
a spacing of 2° . The primary extraction aperture is shown as a solid line and the background annulus as a dashed line. 



confidence level (c.L). iKogut et all |2007h used the full-sky 
WMAP 3-year data to constrain the polarization fraction 
of the AME as traced by FIR maps and concluded that it 
contributes l ess than 1 % of t he ob served polarization sig- 
nal variance. iDickinson et al.1 (^POd ) gave an upper limit of 
10 % (2a) to the possible ex cess emission from th e HII re- 
gion LPH96 at 31 GHz while ICasassus et al.l (|2007l ) gave an 
upper limit of 1 2 % for the excess 31 GHz emission from 
the Helix nebula. M ason et al. | (|2009l ) placed an upper limit 



of 2.7 % (95 % c.L) at 9. 65 GHz for LDN1622 at an angu- 
lar scale of 1.3 arcmin. Ca sassus et al.l (|2008l ) gave a 3<r 
upper limit of 4.8 % at the peak intensity of the p Ophi- 
uchi CB I map at 31 GHz, on s c ales o f « 9 arcmin. More 
recently, lLopez-Caraballo et all (|201ll ) used the WMAP 7- 
year data to place upper limits of 1.0, 1.8 and 2.7%, at 23, 
33 and 41 GHz, respective ly (with a 95% e l. ), wit hin the 
Perseus Molecular Cloud. iMacellari et al.l (|201lh used 
a template fitting technique to constrain the dust- 
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correlated component at 23 GHz to be less than 5 %. 

These measurements therefore indicate that the AME has 
a relatively low (< 3 %) level of polarization, which is con- 
sistent with expectations for spi nning dust models based o n 
resonance relaxation proposed bv lLazarian fc Draind (2000). 
At the same time , they largely rule out m odels of magnetic 
dipole radiation (jDraine fc Lazarianlll999T l. which typically 
predict high polarizatio n fractions (up to ~ 40 %) . However, 
it should be noted that lDraine fc Lazarianl (|l999h also show 
a model with random inclusions of metallic Fe that produces 
very little polarization (< 1%). 

In this Letter, we examine the WMAP 7-year data 
to constrain the polarization fraction for the AME in the 
p Ophiuchi and Perseus clouds. Both regions were stud- 
ied in intensity using data from the Planck satellite and 
were shown t o be dominated by AME at fre quencies of 
~ 20-60 GHz (|Planck Collaboration et alj|201j-h . with spec- 
tra that can be well-fitted by models of spinning dust grains. 
Section [5] describes the WMAP maps used in the analysis. 
Section [3] presents the analysis and a discussion of the re- 
sults. Section |3] concludes. 



2 MAPS 

The analysis is based on the WMAP 7-year data available 
from the LAMBDA websit«fl We used the T-smoothed 
maps available in HEALPix format with N s ide = 512 
i|G6rski et alj|2005h . The maps contain the thermodynamic 
brightness temperature (iiiKcmb) with respect to the CMB 
temperature, at each pixel in Stokes /, Q and U. The data 
also contain estimates of the white noise level in each pixel, 
including the noise correlation between Q and U. 

Fig. [1] shows 10° x 10° WMAP maps of the p Ophi- 
uchi region centred at Galactic coordinates (l,b) — 
(353°.05, +16.°90) at each of the five WMAP bands (see cap- 
tion of Fig. [T]). The maps depict Stokes / (total intensity) 
and linear polarization in the form of Stokes Q and U maps. 
We also show the polarized intensity (P) and polarization 
angle (PA) maps derived directly from the Q and U maps, 
where P = ^/Q 2 + U 2 and PA = § tan _1 ([//Q). In total- 
intensity, there is a bright feature seen in the centre, which 
corresponds to th e p Ophiuchi molecular cloud, or AME- 
G353.05+16.90 of I Planck Collaboration et all l|201lh . It is 
detected with high signal-to-noise ratio at all frequencies, 
with the majority of the flux falling within the 2° aperture. 
There is fainter extended emission around the main feature 
and also in the surrounding region. 

In polarization, the maps look very different. There are 
large-scale features running across the maps, particularly at 
the lower frequencies of K- and Ka-bands. There is a rel- 
atively bright synchrotron spur running diagonally across 
the top left-hand side of the polarization maps, which is 
not obvious from the total-intensity maps. This is presum- 
ably because there is strong AME and free-free emission 
which is dominating over the synchrotron component. How- 
ever, in polarization, the synchrotron component is strong 
while the AME/free-free components is clearly much weaker. 
Away from the bright spur, there is weak polarized emission 

1 http://lambda.gsfc.nasa.gov/ 
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Figure 2. T-T plot of Ka-band (33 GHz) vs K-band (22.7 GHz) 
for a 5° X 5° region centred at (I, b) = (347.°8, +16.°8). The least- 
square straight line is shown. The slope corresponds to a temper- 
ature spectral index = —2.59 ± 0.13. 

in the background of p Ophiuchi region and extending to 
the right-hand side of the map. This produces coherency 
in the polarization angle at K- and Ka-bands (see Fig. [1]). 
A T-T analysis (Fig. [2} of the a 5° x 5° region centred at 
(l,b) = (347.° 8, +16° 8) finds a spectral index between 22.7 
and 33.0 GHz of /3 » -2.6 (T oc v 13 ). This is consistent with 
non-thermal synchrotron emission. At Q-band (40.6 GHz) 
the polarized emission is weaker and the maps are domi- 
nated by noise. 

In addition to the diffuse synchrotron emission, compact 
(polarized) emission may be present in the form of extra- 
galactic radio sources; however, a search of the NASA Extra- 
galactic DatabassQ found no bright (> 100 mjy at WMAP 
frequencies or below) radio sources within a 2° radius that 
would be detected in the WMAP data. 



3 ANALYSIS AND DISCUSSION 

We carried out aperture photo metry of the maps using the 
same 2° diameter aperture as in lPlanck Collaboration et al.l 
(2011), indicated in Fig. [T] The results are given in Ta- 
ble [3] The results for the Perseus cloud, with the same 
aperture centred at (l,b) — (160.°26, — 18.°62), are also 
given. The temperatures were converted to flux density 
units (Jy) for the effective frequencies of WMAP. An es- 
timate of the background was subtracted based on the me- 
dian value of an annulus with an inner radius of 80 arcmin 
and outer radius of 120 arcmin (Fig. [TJ. We also list the 
total-intensity flux density of the spinning dust, I s( j, from 
IPIanck Collaboration et al.l (|201ll ). which appears to domi- 
nate the K-, Ka-, and Q-bands. For Perseus, we could not use 
the flu x densities directly from IPIanck Collaboration et al.l 
l|201lf) because they filtered the data; we re-calculated the 
flux densities for a 2° aperturs[f|. Colour corrections were 
applied to the flux densities for the intensity, using a model 
consisting of free-free, spinning dust, a CMB fluctuation, 
and thermal dust with fixed emissivity index and dust 

2 http://ned.ipac.caltech.edu/ 

3 The flux densities for Perseu s are lower than those of 
IPIanck Co llaboratio n et al. I j201ll) since they used a Gaussian 
model with a larger area than we used here. 
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Table 1. Aperture photometry of the WMAP 7-year maps in the p Ophiuchi (fop) and Perseus (bottom) molecular cloud regions, using 
a 2° diameter aperture. The columns contain flux densities for intensity 7, the spinning dust contribution to the intensity (/ sc j), Q, 
U, the observed polarized intensity (P), the maximum likelihood value for the noise-bias corrected polarized intensity (Po), the total 
polarization fraction (II), and the polarization fraction for spinning dust (n s( j). Upper limits are at the 95 c.l. 



Frequency 
(GHz) 


I 

(Jy) 


(Jy) 


Q 

(Jy) 


U 

(Jy) 


P 

(Jy) 


Po 

(Jy) 


n 

(96) 


n sd 

(%) 


22.7 
33.0 
40.7 
60.6 
93.4 


(K-band) 

'Ka-band) 

[Q-band) 

(V-band) 

(W-band) 


26.3 ±5.5 
30.7 ±5.3 
27.7 ±4.6 
26.3 ±4.5 
63.6 ±8.9 


24.8 ±6.6 
27.2 ±6.3 

21.9 ±5.8 
9.8 ±6.5 

6 ± 15 


0.091 ± 0.096 
-0.27 ±0.12 
-0.07 ±0.18 
0.35 ±0.41 
2.8 ± 1.0 


0.23 ± 0.14 
0.02 ± 0.15 
-0.23 ±0.24 
0.41 ± 0.40 
0.0 ± 1.7 


0.25 ±0.13 
0.27 ±0.12 
0.24 ±0.23 
0.54 ±0.41 
2.8 ± 1.0 


0.21+°;^ (< 0.43) 
0-2it° Q {l (< 0.44) 
0.00+° :30 (< 0.57) 
0.00+°- 61 (< 1.1) 
2.6±{-° 


< 1.6 

< 1.4 

< 2.1 

< 4.2 
4.1 ± 1.8 


< 1.7 

< 1.6 

< 2.6 

< 11 


22.7 
33.0 
40.7 
60.6 
93.4 


(K-band) 

(Ka-band) 

'Q-band) 

(V-band) 

(W-band) 


21.0 ±3.1 
20.4 ±3.0 
16.9 ±2.9 
14.9 ±4.0 
32.4 ±9.8 


16.7 ±3.5 
15.7 ±3.3 
11.6 ±3.3 
5.4 ±4.4 
3 ± 10 


-0.118 ±0.071 
0.02 ±0.13 
0.07 ±0.19 

0.098 ± 0.077 
-0.5 ± 1.3 


0.068 ± 0.069 
-0.06 ±0.15 
-0.25 ±0.21 
-0.64 ±0.59 
-0.3 ± 2.3 


0.136 ±0.070 
0.07 ±0.15 
0.26 ±0.21 
0.65 ±0.39 
0.6 ± 1.4 


0.11llH 6 84« 0-24) 
0.00+ 015 (< 0.30) 
0.00 +0 - 30 (< 0.54) 
0.44±°;£« 1.2) 
0.0+ 15 (< 2.9) 


< 1.1 

< 1.5 

< 3.2 

< 8.1 

< 9.0 


< 1.4 

< 1.9 

< 4.7 

< 22 



temperature derived bv lPlanck Collaboration et alj l|201 ll ). 
These corrections were small (typically ~ 1 % or less). No 
colour corrections were applied in polarization. 

The uncertainty in intensity (7) was estimated using 
the r.m.s. fluctuations within the background annulus, which 
contains a contribution from the instrumental noise and sky 
background fluctuations. The code was validated by calcu- 
lating the flux densities of well-known bright radio sources 
and comparing them to the literature. For Q and U we esti- 
mated the noise-only contribution using monte carlo simula- 
tions. We generated 500 noise-only simulations, based on the 
Q-U noise matrices provided by the WMAP team, smooth- 
ing each noise realisation by the appropriate beam transfer 
function to produce l°-smoothed (Gaussian FWHM) maps. 
In this way, we also take into account the correlated noise 
between the Q and U channels. In the p Ophiuchi region, 
this corresponds to an increase of ~ 30 % in the effective 
noise level. We confirmed that estimates of the background 
fluctuations were similar to those found in the noise-only 
simulations i.e. that we are noise-dominated in polarization. 
We also verified empirically that the background fluctua- 
tions were sub-dominant by repeating the analysis in ten 
random apertures in the vicinity of the p Ophiuchi cloud; the 
scatter in the derived flux densities were consistent (< 3<r) 
with the quoted uncertainties. 

Since we are measuring the polarization properties of 
Galactic signals that are relatively bright in intensity, we 
must be sure that instrumental polarization leakage (from I 
to Q and U), due to bandpass mismatch, is negligible. The 
WMAP I, Q, U maps have been corrected for bandpass mis- 
match by solving for an extra term, S (the "spurious" map), 
which does not vary with parallactic angle (since bandpass 
mismatch only depends on the intensity and spectral shape 
of the signal). To estimate the level of residual leakage we 
measured the polarization of bright Galactic Hn regions, 
using the same technique described above. For M42 (Orion 
nebula), we constrain the polarization leakage to be < 0.1 % 
at K-band and even less in the other channels, and hence 
this effect can be safely neglected. We therefore do not in- 
clude any additional contribution to the uncertainties from 
systematic errors in polarization. 



The observed polarized intensity was calculated as 
P = sjQ 2 + U 2 with an uncertainty propagated from the 
Q and U uncertainties. However, it is well-known that 
P is biased positive due to the polarization noise bias, 
which is particularly importan t at low signal-to-noise ra- 
tios JWardle fc Kronberd Il974l ; ISimmons fc Steward Il985l ; 
IVaillancourtll2006t ). Since it is clear that we are in the low 
signal-to-noise regime, we calculated confidenc e intervals in 
P usi ng a Bayesian technique, as outlined by IVaillancourtl 
(2006), and generalized to allow for non-equal noise values 
in Q and U. The noise estimates from the monte carlo simu- 
lations were used to calculate the maximum likelihood value 
for the polarized intensity (Po), along with 68% and 95% 
confidence interval, correctly taking into account the effect 
of noise bias. The values for both the p Ophiuchi and Perseus 
molecular cloud regions are given in Table [3] The uncertain- 
ties are the 68 % confidence limits around the maximum 
likelihood values while upper limits are given at the 95 % 
c.l. 

From the derived polarization intensity values (Po), 
we can derive the polarization fractions for the total in- 
tensity (II = I/Po) and for the spinning dust component 
(n s d = ^sd/Po); these are listed in Table [3] It can be seen 
that for most channels, the polarization is an upper limit 
at the 95 % c.l. For p Ophiuchi, the upper limits are at the 
< 2 % level at 20-40 GHz where spinning dust is the dom- 
inant component in intensity; the tightest constraint is at 
33 GHz with an upper limit of 1.6 %. At W-band (93.4 GHz) 
there appears to be a hint of a polarization signal at ~ 2.3 a, 
which is visible in the Q map (Fig.[T]). It is not clear whether 
this is a real detection, since W-band is more susceptible 
to 1// noise, and therefore the noise bias/confidence inter- 
val may not be exactly correct. Nevertheless, the Rayleigh- 
Jeans tail of thermal dust is do minant at this frequency 
ijPlanck Collaboration et al.ll201ll ) and we may be detecting 
the polarization of thermal dust emission, where a f» 4 % 
polarization fraction is reasonable. For the Perseus cloud, 
our upper limits of 1.4%, 1.9% and 4.7% at K-, Ka-, and 
Q-bands, respectively, are similar to the values (1.0 %, 1.8% 
and 2.7%) derived bv lLopez-Caraballo et al.l (|201ll ). 

Such a low level of polarization of spinning dust emis- 
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sion is important for two main reasons. Firstly, if this can be 
shown to be true for diffuse emission at high Galactic lati- 
tudes, it will simplify component separation for future ultra- 
sensitive CMB polarization experiments. The dominant po- 
larization foreground at frequencies < 100 GHz will be syn- 
chrotron emission, which is polarize d at the > 20 % level 
at high latitudes (|Kogut et al.l [20071 ). Secondly, the polar- 
ization fraction can rule out some models of dust emission. 
Spinning dust emission is expected to be largely unpolar- 
ized, with predictions of < 1 % at frequencies > 20 GHz 
|Lazarian fe Draind |2000| ) . Magnetic dipole radiation can 
exhibit high (> 10 %) polariz ation at 20-40 GHz when there 
is a single magnetic domain (jDraine fe Lazarianl ll999). This 
is clearly incompatible w i th the constraints derived here; 
see lLopez-Caraballo et alj (|201 lj ') for a discussion. However, 
we note that magnetic dipole radiation may also be rela- 
tively unpolarized (< 1 %) if the magnetic inclusions are 
randomized on the grain surface. Also, the constraints 
on the polarization of AME obtained for molecular 
clouds may not be the same for other interstellar 
phases. For instance, hypothetically, one may have 
single magnetic dipole aligned grains in warm and 
hot interstellar phases, but they coagulate in molec- 
ular clouds to produce large grains with random in- 
clusions. 



4 CONCLUSIONS 

We have searched for polarized radio emission associated 
with the p Ophiuchi molecular cloud using WMAP data. 
The bulk of the emission at 20- 40 GHz is due to AME, most 
likely from spinning dust grains dPlanck Collaboration et al.l 
l201ll ). Thus measurements of the polarization correspond to 
a constraint on the polarization of spinning dust emission. 
We find no significant polarization within an aperture of 2° 
diameter. The upper limits on the fractional polarization 
of spinning dust in the p Ophiuchi cloud are 1.7%, 1.6% 
and 2.6% (at 95% c.l.) at K-, Ka- and Q-bands, respec- 
tively. Such low levels of polarization rule out single domain 
magnetic dipole radiation as the dominant emission mech- 
anism for the 20-40 GHz signal but is consistent with spin- 
ning dust emission. Our results for the Perseus cloud are 
at a similar level and are co mparable to those derived by 
lLopez-Caraballo et ail (|201ll ). If AME at high Galactic lat- 
itudes also exhibits low (< 1 %) levels of polarization then 
the dominant diffuse foreground at frequencies < 100 GHz 
will be synchrotron, thus simplifying component separation 
for CMB polarization studies. 
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